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Abstract 

We discuss the nuclear EMC effect with particular emphasis on recent data for light nuclei 
including ^H, ^He, '^He, ^Be, ^^C and ^^N. In order to verify the consistency of available data, 
we calculate the deviation between different data sets. We find a good agreement between the 
results from the NMC, SLAC E139, and HERMES experiments. However, our analysis indicates 
an overall normalization offset of about 2% in the data from the recent JLab E03-103 experiment 
with respect to previous data for nuclei heavier than "^He. We also discuss the extraction of 
the neutron/proton structure function ratio F2 / F2 from the nuclear ratios ^He/^H and ^H/^H. 
Our analysis shows that the E03-103 data on ^He/^H require a renormalization of about 3% in 
order to be consistent with the F2/F2 ratio obtained from the NMC experiment. After such a 
renormalization, the ^He data from the E03-103 data and HERMES experiments are in a good 
agreement. Finally, we present a detailed comparison between data and model calculations, which 
include a description of the nuclear binding, Fermi motion and off-shell corrections to the structure 
functions of bound proton and neutron, as well as the nuclear pion and shadowing corrections. 
Overall, a good agreement with the available data for all nuclei is obtained. 
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I. INTRODUCTION 

The nuclear EMC effect in the deep-inelastic scattering (DIS) has been discussed since 
early 1980s starting from the observation of a dramatic change in the structure function 
of a heavy nucleus relative to that of the deuteron [1]. The nuclear effects in DIS were 
experimentally measured in the form of the ratio 71{A/B) = F^/F.p of the structure func- 
tions (or the cross sections) of two nuclei (usually a complex nucleus to deuterium) in the 
experiments at CERN [2-6], SLAC [7], DESY [8], FNAL [9]. The data were taken for differ- 
ent nuclear targets and at different regions of the Bjorken x and the invariant momentum 
transfer squared and significant nuclear effects were observed (for a review see [10, 11]). 
Recently, the EOS- 103 experiment at Jlab published a high-statistics measurements of the 
EMC effect for ^He, ^He, ^Be, and ^^C, nuclei emphasizing the region of large x [12]. In this 
article we discuss the statistical consistency of data sets from different experiments collected 
with the same nuclear targets in the region 0.1 < x < 0.7, focusing on the kinematical region 
and the targets used in the E03-103 experiment. We also study the sensitivity of the data 
on the EMC effect in '^He to the modeling of the ratio F^/F^. 

A quantitative understanding the nuclear effects in DIS is important for a number of 
reasons. A proper interpretation of experimental data can provide valuable insights into the 
origin of nuclear force and helps us to understand how the properties of hadrons modify 
in a nuclear medium. It should be also noted that the nuclear data often serve as the 
source of information on hadrons otherwise not directly accessible. A typical example is 
the extraction of the neutron structure function which is usually obtained from deuterium 
and proton data and requires a detailed knowledge of nuclear effects [13]. Other examples 
include the using of charged-lepton and neutrino nuclear data in global QCD fits aiming 
to better determine the proton and neutron parton distribution functions and the higher 
twist terms [14-16]. Understanding the nuclear effects is particularly relevant for precision 
measurements in neutrino physics, where heavy nuclear targets are used in order to collect 
a significant number of interactions. 

A quantitative model for nuclear structure functions was recently developed in Refs.[17, 
18]. This approach accounts for a number of nuclear effects including nuclear Fermi mo- 
tion and binding (FMB), nuclear pion excess, shadowing, and off-shell correction to bound 
nucleon structure functions. The detailed analysis of data published before 1996 on the 
ratios TZ{A/B) for a wide region of nuclear targets shows a good performance of the model 
which was able to describe the observed x,Q^, and A dependencies [17]. In the present 
article we compare the predictions of our model with the recent data from HERMES [8] and 
E03-103 [12] experiments and discuss the role of different nuclear corrections in light nuclei. 

The article is organized as follows. In Sec. II we outline the model used in our studies. In 
Sec. Ill we discuss the data from different experiments and also confront the results of our 
calculation with data. In Sec. IV we summarize the results. 

II. THEORETICAL MODEL 

We recall that while DIS is characterized by a large invariant momentum transfer Q ^ M, 
where M is the nucleon mass, the characteristic longitudinal distance in the target rest 
frame L ~ 2qQ/Q'^ = {Mx)~^ is not small in hadronic scale (see, e.g., Ref.[19]). In nuclei, 
the comparison between L and an average distance between bound nucleons rjvTv gives the 
characteristic regions of the Bjorken variable x, which are goverened by different nuclear 
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effects. As tnn ~ 1 2fm, at large values of x we have L <^ tnn and nuclear DIS can 
be approximated by the incoherent scattering off the bound protons and neutrons (impulse 
approximation, or lA). In this region the major nuclear corrections are due to nuclear binding 
and nucleon momentum distribution (Fermi motion). At small values of x the longitudinal 
distance L becomes large and the lA is difficult to justify. In this region the corrections 
due to scattering off nuclear pions (meson exchange currents), as well as the nuclear (anti- 
) shadowing effect due to coherent multiple interactions of hadronic component of virtual 
intermediate boson with bound nucleons, become important. Summarizing, for the nuclear 
structure function we have different contributions (to be specific we discuss F2, for more 
detail see Ref. [17]): 

F,^ = Fi^^) + 6^F,^ + ^cohF^", (1) 

where the first term in the right-hand side stands for the impulse approximation, and (5,r-p2 
and 5coh-^2 are the corrections due to scattering off the nuclear pion (meson) field and the 
coherent interaction of the intermediate virtual boson with the nuclear target, respectively. 

The lA term dominates at large x. This term can be written as a sum of the proton 
(r = p) and the neutron (r = n) contributions [17]: 

j'F,\x,Q^) =J2f m'P^ie,p) (1 + ^) (y^ + ^) F,^{x',Q\p% (2) 

where the integration is taken over the four- momentum of the bound nucleon p = {M + e, p) 
and [dp] = de d^p/(27r)^. In the integrand, V^^^\e,p) is the proton (neutron) nuclear 
spectral function, which describes the energy and momentum distribution of bound nucleons, 
F^^""^ is the structure function of bound proton (neutron), which depends on the Bjorken 
variable x' = Q'^/{2pq) the momentum transfer square and also on the nucleon invariant 
mass squared = (M + e)^ — p^. In Eq.(2) we use the coordinate system in which the 
momentum transfer q is antiparallel to the z axis, p± is the transverse component of the 
nucleon momentum, and 7^ = 1 + 4x^M^/Q^ and 7'^ = 1 + Ax''^p'^ /Q"^ . 

Nuclei typically have different proton and neutron numbers. For this reason the nuclear 
structure functions generally have both the isoscalar and the isovector contributions. In order 
to separate the isoscalar and isovector contributions in Eq.(2), we introduce the isoscalar 
-pp+n _ -pp j^-pn ^YiQ isovector = ^pp — jy^ spectral functions. The spectral function 
pp+n jg normalized to the total nucleon number A = Z + N with Z and being the proton 
and neutron number, respectively, and P^"" is normalized to Z — N. We separate the 
normalizations from the spectral functions and write 

-pp-n = - N)Vi, (3b) 

where the reduced spectral functions Vq and Vi are both normalized to unity. Using Eq.(3) 
we explicitly write the nuclear structure function in Eq.(2) in terms of the isoscalar and the 
isovector contributions 

F,^/A = (F,^)^ + ^ {FD, , (4) 

where F.^ = |(-F2 +-^2*) is the structure function of the isoscalar nucleon and -Ff "'^ — F2—F2 
and the parameter /3 = {Z — N)/A describes the fractional difference of protons and neutrons 
in a nucleus. The quantities (F)q and {F)-^ are the contracted notations of the integration 
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in Eq.(2) taken with the reduced spectral functions Vq and Vi, respectively. The model of 
Vo and Vi, which is used in this article for nuclei with A > A, is discussed in Sec. II A. We 
also note that for ^He we apply Eq.(2) with the proton and neutron spectral functions taken 
from microscopic calculations of Refs.[21, 23]. 

It should be noted that the experimental data on the nuclear EMC effect have been 
often presented for the isoscalar part of the nuclear structure functions in order to facilitate 
comparison between different nuclei. In order to separate the isoscalar part from nuclear 
data, a common practice is to multiply the data by a factor 



that accounts for the proton- neutron difference in a nuclear target. Indeed, neglecting nu- 
clear effects we have from Eq.(4) CigF.^ = AF.^ . It must be emphasized, however, that this 
method is only an approximate way to isolate the isoscalar contribution to the structure 
functions, as a correct procedure should involve the details of the proton and neutron distri- 
butions in nuclei. Also, it is important to realize that the data are biased by this procedure 
as the factor Cis depends on the neutron structure function through the ratio F2 /F2 and 
different experiments employ different models of this ratio. For this reason the comparison 
between theoretical calculations and data is not clear-cut for nonisoscalar nuclei. Appar- 
ently, a consistent comparison between data and calculations for nonisoscalar nuclei should 
involve unbiased data. However, in practice it is not always possible to remove the isoscalar 
correction from data of specific experiments. 

A. Nuclear Spectral Function 

The nuclear spectral function V describes the energy and momentum distributions of 
bound nucleons and can be written as 



where a^(p) and a{p) are creation and annihilation operators of the nucleon with momentum 
p, Ea is the energy of the nuclear state \A), and H is the Hamiltonian of the system. Note 
that we discuss unpolarized scattering and in Eq.(6) we implicitly assume the sum over 
nucleon polarizations. For simplicity, we also suppress explicit notations for the nucleon 
isospin state (proton or neutron) of V and the operators and a. 

Note also that Eq.(6) is written in the target rest frame and defines the nuclear spectral 
function as a function of nucleon energy e = po — M and momentum p. However, in the 
literature the spectral function is usually considered as a function of the nucleon separation 
energy E. The spectral function in this case is denoted by P{E,p). The relation between 
V{e,p) and P{E,p) is driven by a relation between e and E, which differ by the kinetic 
energy of a recoil nuclear system of A—1 nucleons, e = —E — p^ /2Ma-i (for more details 



A common way to calculate the spectral function is to insert a full set of the intermediate 
states in Eq.(6) and evaluate the sum of the corresponding nuclear matrix elements. In 
the case of the deuteron, the intermediate states reduce to a free proton or neutron, and 
the spectral function is given entirely in terms of the deuteron wave function square and 





V{e,p) = 2n{A\a\p)5{EA -H- e)a{p)\A) 



(6) 



see Ref.[20]). 
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e = Eo — where Eo ^ —2.22 MeV the deuteron binding energy (the exphcit form 

of Eq.(2) in this case is given in Ref.[17]). 

For a ^He nucleus, the proton spectral function receives two contributions: from the 
bound (pn) intermediate state corresponding to a deuteron, where the separation energy 
IS E = Ed — EiYici with Ei^ic ~ —7.72 MeV the ^He binding energy; and from the (pn) 
continuum scattering states. The neutron spectral function, on the other hand, has only the 
{pp) continuum contribution: 

P^{E,p) = r^''\p)6 {E + EsHe - Eo) + F^'^^'\E,p), (7a) 
P"(E,p) = /"(™'^*)(E,p). (7b) 

A number of calculations of ^He spectral function is available. In Ref.[21] this spectral 
function was obtained by solving the Faddeev equation with the Paris NN potential [22] for 
the ground state wave function, and constructing its projection onto the deuteron and two- 
body continuum states. In Ref.[23] the spectral function was obtained using a variational 
^He wave function calculated for a more recent AV14 NN potential. While in our numerical 
analysis we use the spectral functions from Ref. [21], we verified that both Ref. [21] and 
Ref.[23] lead to a consistent 7?.(^He/^II) ratio. ^ 

For A> 4 nuclei we follow the model of Ref. [17], in which the full nuclear spectral function 
was calculated as a sum of the mean-field spectral function, describing a low-energy part 
of the separation energy spectrum, and the term generated by short-range A^A^ correlations 
in the nuclear ground state responsible for a high separation energy and high-momentum 
component. 



Nucleus 


Ea/A 




(p2)/2M 




-1.11 


-11.46 


9.24 


^He 


-2.57 


-17.95 


12.87 


^He 


-7.07 


-40.06 


25.01 


9Be 


-6.46 


-41.20 


27.40 


12c 


-7.68 


-45.35 


28.83 


14n 


-7.48 


-45.13 


28.40 



TABLE I. The nuclear binding energy per nucleon Ea/A, a bound nucleon energy e and kinetic 
energy p^/2M averaged with the nuclear spectral function normalized to one nucleon (all in MeV 
units). 

In order to illustrate the evolution of the nuclear binding effect in light nuclei, in Table I 
we list the values of nuclear binding energy along with the average separation and kinetic 
energies of a bound nucleon for a few light nuclei. For the deutron we use the Paris wave 
function, the parameters for ^He were calculated using the spectral function of Ref. [21], 
whereas for A > 4 nuclei a model spectral function of Ref. [17] was used. Note a dramatic 
change in the energy parameters when going from ^He to ^He, which is also the underlying 
reason of a difference in the magnitude of nuclear corrections to the structure functions, as 
will be discussed in Sec. III. 

^ In fact, the spectral functions of Ref. [21] and Ref. [23] result in an almost identical EMC ratio for x < 0.85. 
For X > 0.85 and values of the E03-103 experiment, the calculation with the spectral function of Ref. [23] 
gives 7^(^He/^H) larger by some 1-2%. For the spectral function of Ref.[23] we used a more recent results 
with AV18 and Urbana IX potential (G. Salme, private communication). 
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B. Nucleon Structure Functions 



In the DIS region, motivated by the twist expansion, the nucleon structure functions can 
be written as 

F,{x,Q') = Fr''{x.Q') + Hi'\x,Q')/Q' (8) 

where F™*^ is the leading-twist (LT) structure function corrected for the target mass effects 
and ifg^^ is the function describing the contribution of the twist-four terms. The LT structure 
functions are computed using the proton and neutron PDFs extracted from analysis of DIS 
and Drell-Yan data in a kinematical region of > 1 GeV^ and invariant mass W > 1.8 GeV 
[24, 25] with the coefficient functions calculated to the NNLO approximation in the strong 
coupling constant in pQCD [26]. The target mass correction is computed following Ref. [27]. 
Although in general the twist expansion should include an infinite chain of the HT power 
terms, recent phenomenology suggests that Eq. (8) with only twist-four correction provides a 
good description of data down to Q ~ 1 GeV [15, 24, 25, 28]. It is also worth noting that this 
model is consistent with the duality principle and on average describes the resonance data 
with W < 1.8 GeV [15], which is relevant for the kinematical region of the JLab E03-103 
experiment [12]. 

In calculating the nuclear structure functions one has to deal with the structure functions 
of the bound proton and neutron, which generally differ from those of the free proton and 
neutron. The effect of modification of bound nucleon SF is related to the dependence on 
the nucleon invariant mass in the off-shell region [see Eq.(2)]. One can separate the two 
sources of dependence: (i) The terms originating from the target mass effect in the 
off-shell region. The corresponding correction is of the order p^/Q"^. In Ref. [17] this effect 
is evaluated by the replacement — )■ p"^ in the expressions of Ref. [27]. (ii) The off-shell 
dependence of the the parton distributions (or LT structure functions); the latter effect is 
generally not suppressed by the inverse powers of Q^. In the description of this effect we 
treat the nucleon virtuality v = {p"^ — M'^)/M'^ as a small parameter, expand in series in v 
keeping the leading term: 

Fi^{x, Q\p') = F^-'ix, Q') (1 + 5Mx, Q') v) , (9) 
6f2 = d\nF^^/dlnp^, (10) 

where the first term on the right in Eq.(9) is the structure function of the on-mass-shell 
nucleon and the derivative is evaluated at p"^ = M^. 

Although the off-shell function (10) generally depends on the structure function type, it 
was suggested in Refs.[17, 18] that the relative off-shell effect (10) is common for all types 
of the nucleon PDFs. Thus, the function Sf2 measures a relative response of the nucleon 
parton distributions to the variation of the nucleon mass and, after the averaging with the 
nuclear spectral function, it describes the modification of the nucleon PDFs in the nuclear 
environment. 

III. COMPARISON WITH DATA 

We start this section by summarizing the results of our previous analysis [17] of data on 
the ratios of structure functions 1Z{A/B) = F^/F^ published before 1996. In Ref. [17] we 
performed a statistical analysis by calculating = ~ T^^^)^ / o'^ {'J^'^'^^) , where TZ^^"^ 
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Targets 






xVdof 




NMC 


EMC 


E139 


E140 


BCDMS 


E665 


HERMES 


^He/^H 


10.8/17 




6.2/21 










^Li/^H 


28.6/17 














9Be/2H 






12.3/21 












14.6/17 




13.0/17 










9Be/i2c 


5.3/15 














i2c/7Li 


41.0/24 














14N/2H 














9.8/12 


27A1/2H 






14.8/21 










27A1/C 


5.7/15 














40Ca/2H 


27.2/16 




14.3/17 










^OCa/^Li 


35.6/24 














40Ca/i2c 


31.8/24 










1.0/5 




56Fe/2H 






18.4/23 


4.5/8 


14.8/10 






56Fe/i2c 


10.3/15 














63Cu/2H 




7.8/10 












84Kr/2H 














4.9/12 


108Ag/2H 






14.9/17 










ii9Sn/i2c 


94.9/161 














197Au/2H 






18.2/21 


2.4/1 








207pb/2H 












5.0/5 




207pb/12c 


6.1/15 










0.2/5 





TABLE II. Values of x^/DOF between different data sets with > 1 GeV^ and the predictions 
of Ref.[17]. The normahzation of each experiment is fixed. The sum over ah data results in 
xVDOF = 466.6/586. 



and TZ are the data points and model predictions, respectively, and a is the uncertainty 
of the experimental data points. The examined data come from different experiments for 
a variety of targets ranging from ^He to 208p]^ with > 1 GeV^ and in a wide region of 
Bjorken x (for a complete list of the analyzed data see Table 1 of Ref. [17]). 

The ratio TZ^^ was calculated using approach outlined in Sec. II. We note that the cal- 
culations with Eq.(2) with the realistic nuclear spectral function, accounting for the Fermi 
motion and nuclear binding effect (FMB), and no off-shell correction explains a general 
trend of the EMC effect at large x values [29-31]. Nevertheless, the FMB effects alone fail 
to quantitatively describe data. In Ref. [17] a hypothesis was tested that the account of the 
off-shell effect in the bound nucleon SF would provide an accurate description of data. The 
function 6f2 was assumed to depend only on x which was then treated phenomenologically as 
a polynomial parametrization. Then the parameters, together with their uncertainties, were 
extracted from statistical analysis of data. This approach lead to an excellent agreement 
with all available data on 71{A/B), reproducing the observed x, Q^, and A dependencies. 
In order to demonstrate the consistency between our model and the data, in Table II we list 
the values of x^/DOF obtained from data on different nuclei and from different experiments. 
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A. Data from nuclear targets with A> A 

More recent data on the nuclear EMC effect are available from HERMES experiment at 
HERA [8] and E03-103 experiment at JLab [12]. The HERMES experiment published the 
data for ^He/^H, ^'^N/^H and ^^Kr/^H ratios obtained with a 27.5 GeV positron beam for 
Bjorken x in the range 0.0125 < x < 0.35 and an average ranging from about 0.5 GeV^ 
in a low x region to about 4 GeV^ at large x. The E03-103 experiment published the 
measurements of ^He/^H, ^He/^H, ^Be/^H, and ^^C/^H ratios with a 5.77 GeV electron 
beam for x values ranging from 0.325 to 0.95 and an average from about 3 GeV^ at 
X ~ 0.3 to 6 GeV^ at high x. In Fig.l we summarize these data together with the former 
data by the NMC and SLAG E139 collaborations. 

1. Consistency of different data sets 

In this section we study the consistency of data from different experiments. From Fig.l we 
can conclude that the slopes of the EMG ratio as a function of x from different experiments 
are in a good agreement. Note that the beam energy of these experiments differs significantly. 
As the beam energy determines the range of possible values for each x-bin, the typical 

values of the NMG experiment are significantly higher than those of the SLAG E139, 
HERMES, and JLab E03-103 experiments. Nevertheless, the data show no systematic 
dependence of 7?. in the valence region 0.25 < x < 0.6. 

Contrary to the slope, the overall normalization factors of different experiments are not 
fully consistent. If we consider the ^^G and ^^N nuclei, which have a similar atomic number, 
we observe a good agreement in the normalizations of the NMG, SLAG E139, and HERMES 
experiments. However, the data points of the JLab E03-103 experiment are somewhat 
shifted above the data from other experiments. A similar trend is also present for the other 
nuclei, as can be seen from Fig.l. It is interesting to observe that for the double ratio 
(9Be/2H)/(^2Q/2jj^ ^j^g normalization of the E03-103 data seems to agree with the SLAG 
E139 and NMG data. This fact points toward a common normalization offset for all ratios, 
possibly related to a common source of systematic uncertainty. It should be noted that the 
overall normalization uncertainty is not shown in Fig.l but is summarized in Table HI. The 
NMG experiment achieved the most precise absolute normalization within 0.4%, followed by 
the SLAG E139 experiment. For this reason we try to use NMG data as a reference and we 
evaluate the consistency of other experiments with NMG. From Fig.l we can see that the 
statistical uncertainties of NMG data become large at x > 0.3 and therefore we can have 
an accurate comparison only in the region 0.1 < a; < 0.3. Following this observation, we 
proceed in two steps. First, we study the consistency of SLAG E139 and HERMES data 
with NMG in this overlap region. Then we compare the E03-103 data with SLAG E139 data 
since those two experiments have a broader overlap region for x > 0.3. 

In order to evaluate the compatibility in the overall normalization factors of different 
experiments we follow a approach. We consider the data shown in Fig.l in the region 
0.1 < a; < 0.7 and we rebin the measurements with a common bin size of 0.1/bin (six bins in 
total from 0.1 to 0.7). Within each larger bin, we calculate the weighted average of the data 
points from each experiments with weights given by l/cr^. For each pair of experiments, we 
then calculate the corresponding in the overlap region assuming a fixed normalization, 
as reported by the experiments. The results are summarized in Tables IV and V. From 
Table IV we can conclude that the three experiments NMG, SLAG E139 and HERMES 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Bjorken x 

FIG. 1. (Color online) The 7?.-ratios of ^^N, ^^C, ^Be, ^He with respect to deuterium and the ratio 
^Be/i^C as measured by the NMC [4-6] (full squares), SLAC E139 [7] (fuh circles), JLab E03- 
103 [12] (open circles) and HERMES [8] (full triangles) experiments. Statistical and systematic 
uncertainties are added in quadrature, while the normalization uncertainty is not shown. For the 
E139 and E03-103 experiments the ratio ^Be/^^C is calculated as a double ratio of data on ^Be/^H 
and 12C/2H. 
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Targets 


NMC 


SLAC E139 


HERMES 


JLab E03-103 


3He/2H 






1.4% 


1.84% 


^He/^H 


0.4% 


2.4% 




1.5% 


9Be/2H 


0.4% 


1.2% 




1.7% 


12C/2H 


0.4% 


1.2% 


1.4% (i^N/^H) 


1.6% 



TABLE III. Overall normalization uncertainties on the ratios TZ from different experiments. 



Experiment 


SLAC E139 NMC 


HERMES 




xVDOF for 12C/2H: 




SLAC E139 


0.7/3 


0.2/2 


NMC 




0.1/3 




xVDOF for ^Be/^H: 




SLAC E139 


0.7/3 






xVDOF for "^He/^H: 




SLAC E139 


2.2/3 





TABLE IV. Compatibility between SLAC E139, NMC and HERMES data in the region 0.1 < x < 
0.5. The data shown in Fig.l have been rebinned with a bin size of 0.1 for all experiments. The 
normalization of each experiment is fixed. The normalization uncertainties are not included in the 
evaluation of x^. 

xVDOF between SLAC E139 and JLab E03-103 





^He/^H 


9Be/2H 


12C/2H 


Total 


Default normalization 


4.6/4 


21.7/4 


16.4 /4 


42.7/12 


Normalization factor 0.98 


1.1/4 


6.4/4 


1.3/4 


8.8/12 



TABLE V. Compatibility between SLAC E139 and JLab E03-103 data in the region 0.3 < x < 0.7. 
The data shown in Figure 1 have been rebinned with a bin size of 0.1 for all experiments. In 
the second row all JLab data have been rescaled by an overall factor 0.98. The normalization 
uncertainties are not included in the evaluation of x^- 



Targets 


xVdof 


Normalization 


Normalization 


xVdof 






factor 


uncertainty 


w/o offset 


^He/^H 


11.9/20 


1.9 ±0.3% 


1.50% 


23.6/21 


9Be/2H 


9.2/20 


2.0+°-^% 


1.70% 


23.0/21 


12C/2H 


5.2/20 




1.60% 


27.0/21 


AU data 


26.3/60 


73.6/63 



TABLE VI. Values of xVDOF between JLab E03-103 data [12] with > 2 GeV^ {x < 0.85) and 
our predictions [17]. The overall normalization of each data set has been kept floating while the 
corresponding normalization uncertainty is not included in the x^- For comparison, the last column 
shows the x^/DOF obtained without any normalization offset but including the normalization 
uncertainty to the evaluation of x^- 



10 



appear to be statistically compatible in the absolute normalization. However, the picture 
differs substantially for the EOS- 103 experiment. The JLab kinematics largely overlaps with 
the SLAC E139 experiment and for this reason we calculate the overall between these two 
measurements. The results shown in Table V signal on inconsistency in the normalization 
of data from these experiments. We found that an overall normalization factor of 0.98lo;oo|, 
applied to the ^He/^H, ^Be/^H, and ^^C/^H data of the E03-103 experiment, makes them 
statistically compatible with the SLAC E139 data, lowering by A^^ = 33.9. Therefore, 
in the following comparisons we will apply an overall normalization factor of 0.98 to the 
measurements "^He/^H, ^Be/^H and ^^C/^H from Ref. [12]. 



2. Comparison with model predictions 

Figure 2 shows a comparison of our calculations with HERMES data [8] and also ^^C/^H 
data from NMC experiment. These data emphasize the region of small x where the major 
nuclear correction is due to nuclear shadowing. We found a good agreement of model 
calculations with data over the whole kinematical region of x and Q^. We note that the 
shadowing effect arise in our model due to medium effects in the propagation of hadronic 
component of the virtual photon in nuclear environment. The magnitude of the shadowing 
effect is driven by the total cross section of scattering of hadronic component off a bound 
nucleon. This quantity strongly depends on in the region of low [17], which is the 
underlying reason for the difference in the shadowing correction for the NMC and HERMES 
data in Fig. 2. We note that characteristic are significantly smaller for the HERMES 
experiment. For example, in the x-bin 0.018 the average is 2.9 GeV^ for NMC and 
about 0.66 GeV^ for HERMES. It should be also noted that pQCD and the twist expansion 
approach employed in Eq.(8) is out of the region of applicability for such low values of Q^. 
In order to evaluate the structure functions and their ratios at low Q^, we apply the spline 
extrapolation of Eq.(8) from the region > 1 GeV^ to — using the requirements 
from the conservation of the electromagnetic current [17]. We also note that a bump in the 
EMC ratio between x = 0.1 and 0.3 is due to an interplay between the nuclear shadowing 
and pion correction effects. 

In Fig.3 we show a detailed comparison of our predictions with the E03-103 data for 
the ratios ^He/^H, ^Be/^H, and ^^C/^H. For each x bin, the ratios of nuclear cross sections 
were calculated at the average value provided by the E03-103 experiment [12]. After 
applying the normalization factor of 0.98, required to match the NMC, E139 and HERMES 
data, we observe a very good agreement between our model and the E03-103 experiment for 
X < 0.85. Table VI lists the results of calculation of between the model predictions and 
the data with and without renormalization. As a cross check, we also list the values of a 
normalization factor minimizing if we leave the normalization floating for each of the three 
ratios ^He/^H, ^Be/^H, and ^^C/^H from the E03-103 experiment. Note that the 2% offset 
appears statistically compatible with the overall normalization uncertainties quoted by the 
E03-103 experiment. Indeed, if we keep fixed the JLab normalization and we rather add the 
corresponding normalization uncertainties into the calculation, we obtain a normalized 

of 1.17. We also remark that allowing the normalizations to float reduces the x^/DOF to 
significantly less than unity that may signal an issue with the quoted uncertainties on the 
data. 

In the region of x > 0.85 our calculations somewhat underpredict the E03-103 data. 
We note that in this region the data of Ref. [12] have < 2 GeV^ and therefore lie 
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FIG. 2. (Color online) Ratios of cross sections for positron scattering off ^ N and ^ Kr nuclear 
targets with respect to a deuterium target. Data from the HERMES experiment [8] (full triangles) 
are compared with predictions of Ref.[17] calculated for the same kinematics (solid lines). Statistical 
and systematic uncertainties are added in quadrature, while the normalization uncertainty is not 
shown. The data from the NMC experiment for ^^C target and the results of Ref.[17] (dash-dotted 
line) are also shown in the top panel to illustrate the effect of the dependence at low x values 
(see text). 
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FIG. 3. (Color online) Data on the TZ ratios of ^^C, ^Be, and ^He with respect to deuterium 
compared with predictions of the model of Ref.[17] for the same kinematics. A common normal- 
ization factor of 0.98 has been applied to all data points of Ref.[12], and statistical and systematic 
uncertainties are added in quadrature. The result of a calculation in impulse approximation with 
no off-shell correction is also shown as dashed-dotted line for comparison. 
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within the resonance region, while in our calculations we use a DIS model for the structure 
functions. Although we verified the consistency of low-Q DIS structure functions with the 
duality principle, a more detailed comparison would require explicit resonance model for the 
structure functions. Furthermore, the region of x > 0.85 is very sensitive to the treatment 
of both the bound nucleon momentum distribution and the target mass correction to the 
nucleon structure functions [17]. In our model, we apply the target mass correction following 
a commonly used method of Ref. [27]. However, it is known that the target mass correction 
of Ref. [27] leads to an incorrect behavior of SF in the limit oi x —> 1 that has a significant 
impact on the calculation of nuclear ratio TZ at very large x values. We also note that in 
the detailed comparison with the data the quasielastic contribution, which is not addressed 
in this article, may not be negligible in this region. 

Finally, Fig. 3 shows the results obtained with and without (5/2 = 0) off-shell correction. 
We recall that the off-shell function 5/2 was determined in the former analysis of data [17]. 
We conclude that the off-shell correction plays a crucial role in understanding the data, 
implying the EMC effect is largely driven by the modification of properties of the bound 
nucleons inside the nuclei. 

B. Data from ^He target and the ratio F^/Fg 

Here we discuss nuclear effects in the ^He target and futher address the problem of 
consistency between data from different experiments. The EMC effect in the ^He nucleus 
was measured by HERMES [8] and recently at JLab [12]. We recall that the data are 
presented in a form corrected for the proton-neutron difference in a nucleus, as discussed 
in Sec. II. The isoscalar correction depends on the ratio F2/F2 and for this reason a direct 
comparison of corrected data is biased by the models used in different experiments. The 
uncorrected data on the EMC ratios are available from Ref. [12] and comparison with those 
data allows us to reduce a bias associated with the isoscalar correction. 

In the former discussion we addressed possible inconsistency in the normalization of the 
EMC ratios from different experiments. Below we discuss the relation between the ratios 
7^(3He/2H) = iFi^y^Ff and ni^R/^R) = ^Ff/F^ and F^/F^, and will argue that the 
comparison of the F2/F2 ratio extracted from ^He/^H data of the E03-103 experiment [12] 
and the ^H/^H data by NMC [32] provides a sensitive test of the normalization of the ^He/^H 



In the absense of nuclear effects, the ratio F2/F2 can directly be calculated from either 
the ^H/^H or '^He/^H ratio. A realistic analysis should include the treatment of nuclear 
effects. Let us first consider the ratio ^H/^H. In order to address the nuclear corrections in 
this problem, it is convenient to consider the ratio R2 = F^ / (F^ + F^) and recast F^ in 
terms of i?2- We have for F^/Ff : 



Addressing the ratio ^He/^H, we need to consider = F2^^/{2F2 + F^) along with R2. 
We recast F2^^ in terms of R^ and after a simple algebra we have: 



data. 



F^/F^ = 27^(2H/^H)/i^2 - 1. 



(11) 




(2-z)/(.-l) 



(12) 



where z 



|7^(3He/2H)/^2/i^3• 
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In Fig. 4 we show the results of extraction of jFl^ from the NMC and the E03-103 data 
using Eqs. (11) and (12). It should be remarked that a typical region of is significantly 
higher in the NMC experiment than that of the E03-103 experiment at JLab. Neverthe- 
less, Ref. [32] presents the results for the dependence of the ^H/^H ratio. In order to 
reduce a bias due to dependence, in our analysis we use the NMC parametrization of 
dependence of their data and select the values of close to those of Ref. [12]. 

In order to test the impact of nuclear effects on the extraction of jFl^^^ we performed the 
analysis with the full treatment of nuclear effects and also with no nuclear effects included, 
i.e., assuming R2 = = 1. Both results are shown in Fig. 4. The ratios R2 and i?3, used 
in the full analysis, are shown in Fig. 5. The ratio R2 was calculated as described in Ref. [17] 
using the proton and neutron structure functions of Refs.[15, 33] and the Paris deuteron 
wave function [22], while the ^He structure function and the ratio -R3 was calculated by 
Eq.(2) using the spectral function of Ref. [21] (see Sec. II A). We observe from Fig. 5 that 
the EMC effect in R^ is similar to that in R2 with the minimum at x ~ 0.7 being somewhat 
deeper for R^. Note that the nuclear effects cancel in the region x ~ 0.35, which is consistent 
with the measurements of the EMC effect on other nuclei. 

A significant mismatch in the values of F^jFl^ extracted from different experiments is 
indicated in Fig. 4. In particular, at a; = 0.35 the central point of the F^/-F|' ratio from 
^He/^H data needs to be rescaled by a factor O.85I0 03 with respect to the corresponding 
value from the ^H/^H data. In this region the extraction of F2/F2 from nuclear data is 
practically unaffected by nuclear effects, as illustrated in Fig. 4. Note, however, that for 
larger x the nuclear corrections are significant. Curiously enough, the sign of nuclear effect 
on the ratio F2 / F^ extracted from ^He/^H differs from that extracted from the ^H/^H ratio. 
This is because the dependence on R2 in Eq.(ll) and the dependence on R2/R3 in Eq.(12) 
essencially differ. Note also that in the case of ^He/^H the nuclear effects come through the 
ratio R2/R3 thus reducing the impact of the Fermi motion at large x. 

The normalization of F2 /F2 is directly related to the normalization of the '^He/^H data 
and the observed mismatch may signal an inconsistency in the normalization of nuclear data. 
We apply Eq.(12) to evaluate necessary renormalization of 7^(^He/^H) in order to match 
the F2 / F2 ratio from ^H/^H data and find that the ^He/^H ratio should be increased by a 
factor 1.03^[J;J^[!^. 

In order to further study the nuclear data normalization issue, in Fig. 6 we compare the 
^He/^H data from the HERMES [8] and JLab E03-103 [12] experiments. The HERMES data 
have been corrected for the proton excess using a correction factor in Eq.(5) with Fg'" from 
the NMC measurements. The open circles correspond to the "isoscalar" data as provided 
in Ref. [12]. The filled circles stand for the E03-103 data corrected for nonisoscalarity in a 
different way: We first obtain a raw data by removing the isoscalarity correction of Ref. [12] 
and then apply the factor Cis by Eq.(5) with F| and of Ref. [15, 33] calculated for 
each X bin. In addition we also apply a renormalization factor of 1.03 to each data point. 
Also shown are the results of our calculations of Cis7^(^He/^H). It is worth noting that 
the isoscalarity correction factor C^s differs significantly from the corresponding correction 
in Ref. [12] and we illustrate its model dependence in Fig. 7. In order to avoid confusion 
caused by different correction factors, in Fig. 6 we compare the calculations with the data 
corrected in a similar way. We find from Fig. 6 that at the overlap region at x = 0.35, such 
renormalized E03-103 data are in a good agreement with HERMES data. 

Finally, we note that the nuclear effects in ^He are significantly smaller than in ^He (cf. 
Fig.3 and Fig.6). Such a difference is naturally explained in our approach by a dramatic 
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FIG. 4. (Color online) The results of the extraction of the ratio from data on the ^H/^H 

cross-section ratio by the NMC experiment [32] (green squares) and the data on the ^He/^H cross- 
section ratio by JLab E03-103 experiment [12] (red circles). The full symbols are obtained with 
the model of nuclear effects of Ref. [17], while the open symbols show the points obtained with 
no nuclear corrections. The solid curve is the result of a calculation according to Eq.(8) with the 
NNLO PDFs and the HT corrections of Refs. [15, 25, 33], for the corresponding to the NMC 
kinematics. The dashed line is a independent model used in Ref. [7]. 
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FIG. 5. (Color online) The ratios R2 and i?3 calculated at the values of x and of Ref. [12] for 
X > 0.3 and at fixed = 3 QeV^ for x < 0.3. 
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FIG. 6. (Color online) Ratio of cross sections for electron scattering off '^He with respect to 
a deuterium target, corrected for the proton excess. Data from the HERMES experiment [8] 
(triangles) and the JLab E03-103 experiment [12] (open circles) are compared with the E03-103 
raw data renormalized by a factor 1.03 and corrected for the proton excess according to Eq.(5) and 
using the proton and neutron structure functions of Refs. [15, 33] (full circles). Our predictions 
calculated for the same kinematics with the complete model described in this article (solid line) 
and in the impulse approximation (dashed line) are also shown for comparison. See text for details. 
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FIG. 7. (Color online) The isoscalar correction for ^He calculated by Eq.(5) using different models 
ofF^VFf. 
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difference in nuclear binding for these two nuclei, as illustrated in Table I. The off-shell effect 
is also smaller in ^He because its rate is controlled by the nucleon virtuality v [see Eq.(9)], 
whose value is driven by nuclear binding. Using Table I we obtain for average virtuality 
V = —0.065 (^He) and —0.14 (^He). Thus, a significantly larger value of \v\ for ^He explains 
more pronounced off-shell effect in this nucleus. 

IV. SUMMARY 

We presented a detailed analysis of nuclear EMC effect for light nuclei. Focusing on 
the overlap region 0.1 < x < 0.7, we perform a statistical analysis of the consistency 
between the data sets from different experiments. We found a good agreement between 
the normalizations of the NMC, SLAC E139 and HERMES experiments for nuclei with 
A > A. However, the points from the JLab E013-103 experiment appear systematically 
shifted above the SLAC E139 data by an overall normalization factor of 0.981q'qo3, common 
to all discussed nuclei with A > 4. This renormalization factor is statistically consistent 
with the normalization uncertainty quoted in Ref. [12]. After applying the normalization 
offset the JLab E013-103 data are in good agreement with both SLAC E139 and NMC data. 

The predictions of Ref. [17] are in a good agreement with renormalized E013-103 data 
for all studied nuclei, as discussed in Sec.IIIA2. We would like to emphasize an important 
role of off-shell correction, which is crucial in the description of both the slope of TZ{x) for 
0.3 < a; < 0.7 and position of its minimum at a; ~ 0.75. This correction is controlled by a 
product 6f2{x)v averaged with the nuclear spectral function. The function Sf2{x) measures 
a relative change in the nucleon structure function due to variation of its invariant mass. 
In Ref. [17] this function was phenomenologically derived from a data analysis on nuclear 
ratios TZ. Further interpretation of observed behavior of 6f2{x) in terms of detailed models 
would help to better understand mechanisms of modification of the partonic structure of the 
nucleon in nuclear environment. 

The HERMES data on ^^N are consistent with the NMC data on g^t x > 0.1. The pre- 
dictions by Ref. [17] agree with both experiments. We observe that at small x the shadowing 
effect is more pronounced in the HERMES data than in the NMC data. This difference can 
be attributed to the dependence, since the average of the HERMES experiment is 
significantly lower than the corresponding one of the NMC experiment. This effect is also 
confirmed by calculations in a model of Ref. [17]. 

A significant part of the present analysis is devoted to the study of the EMC effect in ^He, 
for which the data from both HERMES and JLab E03-103 experiments are available. In the 
overlap region at x ~ 0.35 the normalizations of the two experiments somewhat disagree. 
It is important to note that the data on nuclear ratios are usually corrected for the proton 
excess, which depends on the ratio FJ^/F^ . For this reason, the data appear to be biased by 
different models of the isoscalarity correction. In order to verify the consistency of data, we 
study the relation between the ratios of nuclear structure functions and F^/Ff . We extract 
F^/F^ from both the E03-103 data on the ^He/^H ratio and the NMC data on the ^H/^H 
ratio. A significant difference between the results of these two extractions is observed. In 
particular, we find that at x = 0.35 and ~ 3 GeV^ the ratio FJ'/F|' obtained from the 
JLab E03-103 data is about 15% larger than that extracted from the NMC data. Both 
extractions of F2 / F^ become consistent once a normalization factor of 1.03^0 gpg is applied 
to the ^He/^H ratio of the E03-103 experiment. After such renormalization the E03-103 and 
HERMES data on the ^He/^H ratio also become consistent, and our predictions are in good 



18 



agreement with both data sets in the region x < 0.85. 
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